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Abstract. We consider the quadratic family of maps given by f a {x) = 1 — ax 2 with x E 
[—1, 1], where a is a Benedicks-Carleson parameter. For each of these chaotic dynamical 
systems we study the extreme value distribution of the stationary stochastic processes 
Xq,Xx, . . ., given by X n = for every integer n > 0, where each random variable 
X n is distributed according to the unique absolutely continuous, invariant probability of 
f a . Using techniques developed by Benedicks and Carleson, we show that the limiting 
distribution of M n = max{A"o, . . . , is the same as that which would apply if the 

sequence Xq, X\, . . . was independent and identically distributed. This result allows us to 
conclude that the asymptotic distribution of M n is of Type III (Wcibull). 



1. Introduction 

The dynamical systems we consider in this work are the quadratic maps given by f a (x) = 
1 — ax 2 on I = [— 1, 1], with a G BC, where BC is the Benedicks-Carleson parameter set 
introduced in [BC85j . The set BC has positive Lebesgue measure and is built in such a 
way that, for every a G BC, the Collet- Eckmann condition holds: there is exponential 
growth of the derivative of f a along the critical orbit, i.e., there is c > such that 

|(/a n )'(/a(0))|>e- 

for all n G N. This property guarantees not only the non-existence of an attracting 
periodic orbit but also the existence of an ergodic / a -invariant probability measure \x a 
that is absolutely continuous with respect to Lebesgue measure on [—1, 1]. In fact, these 
Benedicks-Carleson systems are chaotic and highly sensitive on initial conditions. Actually, 
after some iterates, the behavior of most orbits becomes erratic and distributed on the set 
[—1, 1] according to the invariant measure fi a . Hence, it is meaningful to study the statistical 
properties of the orbits of these systems and, here, we are particularly concerned with their 
extreme type behavior. 

A natural way to build a stationary stochastic process associated to f a , for some a G BC, 
is to consider the random variable (r.v.) X defined on the probability space ([— 1, l],/x a ), 
taking values on [—1,1] with distribution function (d.f.) G a {x) = /i a {(— oo,x] fl [—1,1]} 
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and then iterate it by f a , i.e., define 

X n = f a (X n _ 1 ) = f:(X ), for each n G N. (1.1) 

This way, we obtain a stationary stochastic process X , X\, X 2 , . . . with common marginal 
d.f. given by G a (x) = fi a {Xq < x}. The stationarity results from the / a -invariance of the 
probability measure \i a (see for example |KT66l Section 15.4]). Our goal is to study the 
asymptotic distribution of the partial maximum 

M n = mzx{X ,X 1 ,...,X n _ 1 }, (1.2) 

when properly normalized. 

The study of distributional properties of the higher order statistics of a sample, like 
the maximum of the sample is the purpose of Extreme Value Theory. Its major classical 
theorem asserts that there are only three types of non-degenerate asymptotic distributions 
for the maximum of an independent and identically distributed (i.i.d.) sample under linear 
normalization. 

The main result of this work states that the limiting law of M n is the same as if 
Xq,Xi,... were independent with the same d.f. G a . In fact, we verify that under ap- 
propriate normalization the asymptotic distribution of M n is of type III (Weibull). As 
usual, we denote by G~ l the generalized inverse of the d.f. G a , which is to say that 
G-\y) := inf{x : G a (x) > y}. 

Theorem A. For each a G EC and every stationary stochastic process (Aj) igNo given by 
( II. ip . we have: 

P{a n (M n - 1) < x} -> H(x) = 

where a n = (l - G~ l (l - £)) \ 

Haiman |Hai03j has obtained a similar asymptotic result for the natural stochastic pro- 
cess associated with the tent map. The arguments used here for quadratic maps would 
allow us to obtain a different proof of Haiman's Theorem. 

Also, a study concerning extremes for dynamical systems, essentially focusing the finite 
sample behavior of maxima, has already been done by Balakrishnan, Nicolis and Nicolis 
in [BNN95] . 

More recently, Collet jCoOl] has studied the asymptotic distribution of the partial maxi- 
mum of the stochastic process Yi, I2, • • ■ given by Y n = — log |£ — f n \, where £ is a reference 
point in the phase space and / a C 2 non- uniformly hyperbolic interval map. He obtained 
that the Gumbel's limiting law applies for Lebesgue-a.e. point £ chosen in the phase space. 
However, the critical point and its orbit are not included in this full Lebesgue measure set 
of points. 

2. Motivation and Strategy 

The study of the limit behavior for maxima of a stationary process can be reduced, under 
adequate conditions on the dependence structure, to the Classical Extreme Value Theory 



e"(- x ) 1/2 , x < 
, x > ' 
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for sequences of i.i.d. random variables. Hence, to the stationary process X@^ Xi^ • • • we 
associate an independent sequence of r.v. denoted by Zq, Z±, . . . with common d.f. given 
by G a (x) = P {Xq < x} = n a {X < x}. We also set for each n G N 

M n = wax{Z ,...,Z r ^ l }. (2.1) 

Let us focus on the conditions that allow us to relate the asymptotic distribution of M n 
with that of M n . Following [LLR83j we refer to these conditions as D(u n ) and D'(u n ), 
where u n is a suitable sequence of thresholds converging to max^^^i] X (z) — 1, as n goes 
to oo, that will be defined below. D(u n ) imposes a certain type of distributional mixing 
property. Essentially, it says that the dependence between some special type of events 
fades away as they become more and more apart in the time line. D'(u n ) restricts the 
appearance of clusters, that is, it makes the occurrence of consecutive 'exceedances' of the 
level u n an unlikely event. 

Consider a sequence of stationary r.v. Y 1 ,Y 2 ... with common d.f. F. We say that 
an exceedance of the level u n occurs at time i if Yj > u n . The probability of such an 
exceedance is 1 — F{u n ) and so the mean value of the number of exceedances occurring up 
to n is n(l — F(u n )). The sequences of levels u n we consider are such that n(l — F(u n )) — > r 
as n — > oo, for some r > 0, which means that, in a time period of length n, the expected 
number of exceedances is approximately r. 

Observe that, in the case of our stationary process Xq,X\, . . ., we have Xj > u n if and 

only if < \J l ~^ n ■ Since u n — > 1 as n — > oo, then, for large n, an exceedance of the 

level u n is always preceded by a return to a tight vicinity of the critical point. Hence, as 
one would expect, there is an intimate relation between exceedances and deep returns to 
the vicinity of the critical point. 

D{u n ) is a type of mixing requirement specially adapted to Extreme Value Theory. In 
this context, the events of interest are those of the form {Xi < u} and their intersections. 
Observe that {M n < u} is just {X Q < u, . . . ,X n _i < u}. A natural mixing condition in 
this context is the following. Let ^(xi, . . . , x n ) denote the joint d.f. of X il , . . . , X in 
and set G iu ,„ An (u) = G ilr __ >in (u, ...,u). 

Condition (D(u n )). We say that D{u n ) holds for the sequence A / o,A A i,X2, ... if for any 
integers %\ < . . . < i p and j\ < . . . < jk for which j\ — i p > m, and any large nGff, 

\G il ,...,i p j 1 ,...j i (u n ) - G il ^ iv {u n )G ju ..., Jk {u n )\ < 7 (m), 

where 7(771) — > as m — > 00. 

We remark that the actual definition of D(u n ) appearing in [LLR831 Section 3.2] is a 
weaker requirement. Following |Co01j . we proposed, very recently, in |FF08j . a reformula- 
tion of condition D(u n ) that follows almost immediately from the decay of correlations of 
the system and together with D'{un) still allows to relate the asymptotic distributions of 
M n and M n . 
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Condition (D'(u n )). We say that D'(u n ) holds for the sequence X , Xi, X 2 , ... if 



The sequence u n is such that the average number of exceedances in the time interval 
{0, . . . ,[n/k}} is approximately r/k, which goes to zero as k — > oo. However, the ex- 
ceedances may have a tendency to be concentrated in the time period following the first 
exceedance at time 0. Condition 12.21 prevents this from happening, i.e., forbids the con- 
centration of exceedances by bounding the probability of more than one exceedance in the 
time interval {0, . . . , [n/fc]}. This guarantees that the exceedances should appear scattered 
through the time period {0, . . . , n — 1}. 

The special relevance of both these conditions is connected to Theorem 3.5.2 of [LLR83] : 
let a n and b n be sequences such that P{a n (M n — b n ) < x} — > H(x) for some non-degenerate 
d.f. H; if D(u n ), D'{u n ) hold for the stationary sequence X\, X 2 , . . ., when u n = x/a n + b n 
for each x, then P{a n (M n — b n ) < x} — > H(x). This means that if we are able to show 
that conditions D{u n ) and D'{u n ) are valid for the stationary process X , X±, . . ., then 
M n and M n share the same asymptotic distribution with the same normalizing sequences. 
Consequently, our strategy to prove Theorem |2] is the following: 

• Compute the limiting distribution of M n and the associated normalizing sequences 



• Show that conditions D(u n ) and D'(u n ) are valid for the stochastic process 
X , X x , X 2 ,... defined in (Oil . 

The rest of the paper is dedicated to the proof of these assertions and is structured as 
follows. In Sectional we describe the properties of the dynamical systems f a , with a G BC, 
and the Benedicks- Carleson techniques. The validity of condition D(u n ) is a consequence 
of the very good mixing properties of the systems considered here. Actually, it follows 
from the fact that these systems possess a weak-Bernoulli generator (see Section 13.81 and 
Remark 13. ip . Then, in Section HI we study the asymptotic behavior of the maximum in 
the i.i.d. case, identify the limiting distribution of M n and the respective normalizing 
sequences a n and b n . Hence, we are left with the burden of proving D'(u n ). In Section [51 
we use the geometric properties of the systems to show Proposition 15.21 that paves the way 
for the proof of D'(u n ) which is finally established in Section [61 In Section [3, we present 
a small simulation study in order to compare the finite sample behavior of the normalized 
M n with the asymptotic one. 



The Benedicks- Carleson Theorem (see |BC85j or Section 2 of |BC91j ) states that there 
exists a positive Lebesgue measure set of parameters, BC, verifying 



[n/k] 




(2.2) 



a„ and h 



3. Properties of the Benedicks-Carleson parameters 



there is c> (c w log 2) such that \Df™(f a (0))\ > e cn for all n > 0; 
there is a small a > such that |/„ (0)| > e' aVE for all n > 1. 



(EG) 
(BA) 
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Before we describe the Benedicks- Car leson strategy let us have an overview of its key 
ingredients. The critical region is the interval (—5,5), where 5 = e _A > is chosen small 
but much larger than 2 — a. This region is partitioned into the intervals 

(-5,5) = |J I m , 

m>A 

where I m = (e~( m+1 \e~ m ] for m > and I m = —I~ m for m < 0; then each I m is further 
subdivided into m 2 intervals {I m ,j} of equal length inducing the partition Vo of [—1, 1] into 

[-l,-5)u\Jl miJ U(5,l]. (3.1) 

Given J G V, let nJ denote the interval n times the length of J centered at J and define 
U m '■= (—e~ m ,e~ m ), for every m G N. In order to study the growth of Df£(x), for 
x G [—1,1] and a G BC , the orbit is splitted in free periods and bound periods. During 
the former we are certain that the orbit never visits the critical region. The latter begin 
when the orbit returns to the critical region and initiates a bound to the critical point, 
accompanying its early iterates. The behavior of the derivative during these periods is 
detailed in Subsections 13.11 and 13.21 

3.1. Expansion outside the critical region. There is Cq > and M G N such that 

(1) If x, ... , ft\x) £ (-5, 5) and k > M , then \Df k (x)\ > e c ° k ; 

(2) If x, ... , ft\x) i (-5, 5) and f k a (x) G (-5, 5), then \Df k a (x) \ > e c ° fc ; 

(3) If x, ... , ft\x) i (-5, 5), then \Df k a (x)\ > 5e c ° fc . 

If we were capable of keeping the orbit of x away from the critical region then it would 
be in free period forever and the estimates above would apply. However, it is inevitable 
that almost every x G [—1, 1] makes a return to the critical region. We say that n G N is 
a return time of the orbit of x if f™(x) G (—5, 5). Every free period of x ends with a free 
return to the critical region. We say that the return had a depth of m G N if f™(x) G I± m , 
which is equivalent to saying that m = [— log |J" (x) {]. Once in the critical region, the orbit 
of x initiates a binding with the critical point. 

3.2. Bound period definition and properties. Let (3 = 14a. For x G (—5, 5) define 
p(x) to be the largest integer p such that 

\f k a (x) ~ f k a(0)\ < e-?\ Vk<p. (3.2) 

Then 

(1) \\m\ < p(x) < 3\m\, for each x G I m ; 

(2) \DfP(x)\ > e c ' p , where d = ^ > 0. 

The orbit of x is said to be bound to the critical point during the period < k < p. We 
may assume that p is constant on each I m ,j- Note that during the bound period the orbit 
of x may return to the critical region. These instants are called bound return times. 

Roughly speaking, the idea behind the proof of Benedicks-Carleson Theorem is that 
while the orbit of the critical point is outside the critical region we have expansion (see 
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Subsection 13. 1 p : when it returns we have a serious setback in the expansion but then, 
by continuity, the orbit repeats its early history regaining expansion on account of (1EGI) . 
To arrange for the exponential growth of the derivative along the critical orbit f lEGj) one 
has to guarantee that the losses at the returns are not too drastic; hence, by parameter 
elimination, the basic assumption condition (IBAp is imposed. The argument is mounted in 
a very intricate induction scheme that guarantees both the conditions for the parameters 
that survive the exclusions. The condition (lEGfl is usually known as the Collet-Eckmann 
condition and it was introduced in |GE83j . 

3.3. Bookkeeping, essential and inessential returns. A sequence of partitions Vq -< 
V\ -< . . . is built so that points in the same element of the partition V n have the same 
history up to time n. For a detailed description of the construction of this sequence of 
partitions in the phase space setting we refer to [Fre05t Section 4]. Here, we highlight some 
of the main aspects of its construction. 

For Lebesgue almost every x G I, {x} = n n >oU n (x), where co n (x) is the element of V n 
containing x. For such x there is a sequence t±,t2, ■ ■ • corresponding to the instants when 
the orbit of x experiences a free essential return situation, which means C (u ti -i(x)) 
for some \m\ > A and 1 < k < m? . We have that u n (x) = ut i _ 1 (x), for every < n < U 
and f^{uj ti {x)) = ^(/^(x)), except for the points at the two ends of (oj ti l (x)) for which 
it may occur an adjoining to the neighboring interval. If ij is an essential return situation 
for x, then it is either an essential return time for x, which means that there exists m > A 
and 1 < k < m 2 such that C fai^uix)) C 3I m .k] or an escaping time for x, which 
is to say that /(a-i),i C ft(wu{ x )) c M or J -(A-i),i C ft(wu{ x )) c where 
J±(A-i),i is the subinterval of i±(A-i) closest to 0. 

We remark that every point in u G V n has the same history up to n, in the sense that 
they have the same free periods, return to the critical region simultaneously, with the same 
depth and their bound periods expire at the same time. 

We say that v is a free return time for x of inessential type if f%(u v (x)) C 3I m> k, for 
some \m\ > A and 1 < k < m 2 , but f^(u v (x)) is not large enough to contain an interval 
I mt k for some \m\ > A and 1 < k < m 2 . 

3.4. Distortion of the derivative. The sequence of partitions described above is de- 
signed so that we have bounded distortion in each element of the partition V n -\ up to time 
n. To be more precise, consider u G V n -\- There exists a constant C independent of u, n 
and the parameter a such that for every x, y G u>, 

\Df2(x)\ 



\Df£(y)\ 

See |Fre05( Lemma 4.2] for a proof. 



< C. (3.3) 



3.5. Growth of returning and escaping components. Let t be an essential return 
time for ui G Vt, with C fl{oo) C 3I m ,k f or some m > A and 1 < k < m 2 . If n is the 
next free return situation for u (either essential or inessential) then 

|/;(u)|>e CM e- 5/)|m| , (3.4) 
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where q = n — (t + p). See [Fre05t Lemma 4.1]. 

Suppose that oj G Vt is an escape component. Then, in the next return situation for u, 
at time t\, we have that 

|/^M|>e-^. (3.5) 

See [MS93j . [Fre06| Lemma 4.2] or |Mor93| Lemma 5.1] for a proof of a similar statement 
in the space of parameters. 

3.6. Existence of absolutely continuous invariant measures. For every a G BC, the 
quadratic map f a has an invariant probability measure p a that is absolutely continuous with 
respect to Lebesgue measure on [—1, 1]. The existence of absolutely continuous invariant 
measures (a.c.i.m) for a positive Lebesgue measure set of parameters was first proved 
by Jakobson in |Jak81] and others followed. See, for example, |CE83j . |BC85j . |Now85] . 
|Ryc88| , |BY92] and |You92] . 

The a.c.i.m. fi a = p a dx has the following properties: 

(1) p a is the only a.c.i.m. of f a ; 

(2) (f a , Ha) is exact; 



(3) p a = P1+P2, where pi has bounded variation and < p^x) < const 



(1-9)- 



(4) The support of p a is [fa (0), / a (0)] and inf x . 6[/ 2 (0)Ja(0)] p a (x) > 0. 
The proof of these statements can be found in |You92[ Theorems 1 and 2] . 

3.7. Decay of correlations and Central Limit Theorem. The Benedicks- Car leson 
quadratic maps have good statistical behavior. In fact, L. S. Young proved that these maps 
have exponential decay of correlations and satisfy the Central Limit Theorem ( |You92[ 
Theorems 3 and 4]). This was also obtained by Keller and Nowicki in |KN 92j. To be more 
precise, for every a G BC, there exists q G (0, 1) such that for all (p,if) : [—1, 1] — > R with 
bounded variation, there is C = C(<p, tp) such that 

if-i^O f2)dp a ~ J Vdpa J i)dp a < CV\ Mn > 0. (3.6) 

Moreover, if J (fdp a = then for every x G K we have 

A*.(iX>o£<xj >Q(x/a), (3.7) 

I « /n ' ' n— inn 



i=0 



where we are assuming that a := lim^oo ^= / (Ya=o V ° fa) dp a 
notes the N(Q, 1) d.f. 



1/2 

> and $(•) de- 



3.8. Exponential weak-Bernoulli mixing. Keller [Kel94] has obtained a result even 
sharper than (13.61) . Consider the partition of [—1, 1] given by Q = {[—1,0), [0, 1]}. Also, 
for integers k < I, denote by Q l k the join of partitions \J l i=k f~ l Q and by T l k the cx-algebra 
generated by Q l k . According to [Kel94] the partition Q is a weak-Bernoulli generator for 
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every f a with a G EC. This means that the a-algebra Tq 00 coincides, up to sets of Lebesgue 
measure 0, with the Borel cx-algebra of sets in [—1, 1] and that 

Pm(fa, Q, AO 0, as 771 ->■ 00, 



where 



m {fa, Q,fi a ) ■ = 2 sup / sup {\fi a (A\^) - fi a (A)\ : A G T™ +m } dfx a 

k>0 J 



= sup J2 \ti a {AnB)-fi a (A)fi a (B)\. 

k>l,L>l , — ' 

In fact, |Kel94l Theorem 1] states that there are constants C > and < r < 1 such that 

P m (fa,Q,Va)<Cr m (3.8) 

for all m G N. 

Remark 3.1. We observe that if we refine the partition Q by adding one point so that 
{X > u} G To, where To is the cr-algebra generated by Q, then Keller's argument still 
holds with the same type of estimate as in f)3.8p . As a consequence, condition D{u n ) is 
true for every considered sequence u n . 

4. Domain of attraction of the associated i.i.d. process 

We recall that to every stationary stochastic process X , Xi, X 2 , . . . defined in f 1 1.1 1) we 
associated an i.i.d. sequence of r.v. Z , Zi, Z 2 , . . . with common d.f. given by G a (x) = 
P{Xo < x} = /i a {(— oo,x] fl [—1,1]} (see Section [2]). In this Section we will determine 
the domain of attraction corresponding to the d.f. G a , i.e., we will compute the limiting 
distribution of M n , defined in (12.1ft . when properly normalized. For that purpose one must 
look at the tail behavior of 1 — G a (x) as x gets closer to sup^gjj {G a (y) < 1} = 1. According 
to Section I3TB1 (J3J), if 2; is close to 1 we may write fi a {(z, 1]} x y/1 — z), in the sense that 
M °/T^iP "~ * f° r some c > 0, as z — > 1. Hence, for s > sufficiently close to we have: 



1 - G a (1 - s) x Vl - (1 - 8) X yfi, (4.1) 

which means that lim s ^. + 1 ~ G ^ 1 ~ S ^ > 0. At this point, we apply |LLR83t Theorem 1.6.2] 
to obtain that G a , in this case, belongs to the domain of attraction of type III (Weibull) 
with parameter 1/2, since for every x > 



l-G a (l-xh) Vxh 1/2 

hm — — — — = hm — — = x ' . 

h^o+ 1 — G a (1 — h) h^o+ y/h 

Moreover, according to |LLR83| Corollary 1.6.3], if we consider the sequences defined for 
each n G N by b n = 1 and a n — (1 — G~ 1 (l — 1/n)) , where G' a 1 {y) = inf{x : G a (x) > y}, 
then 

f - 1 f e~ (-a:)1/2 x < 

p[a n (M n -b n )<xj ^H(x) = 



x > 



as n — > 00. 
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5. Probability of an essential return reaching a certain depth 

In the study of extremes, one is mostly interested in the probability of occurrence of 
exceedances of the level u n . As we have already mentioned in Section [2], these events are 
related with the occurrence of deep returns. Thus, in this section we do some preparatory 
work by estimating the probability of the returns hitting a given depth. 

For each x £ /, let v n (x) denote the number of essential return situations of x between 
1 and n, s n (x) be the number of those which are actual essential return times and & n the 
number of the latter that correspond to deep essential returns of the orbit of x, i.e, with 
return depths above a threshold > A. Observe that v n {x) — s n {x) is the exact number 
of escaping situations of the orbit of x, up to n. 

Given the integers 0<s<^,s<t><n and the s integers jx, . . . ,j s , each greater 
than or equal to 0, we define the event: 



yi> (n) = < x £ I : v n (x) = v, & n (x) = s and the depth of the i-th deep essen- > . 
[ tial return is ji Vi £ {1, . . . , s} J 

Remark 5.1. Observe that the upper bound ^ for the number of deep essential returns up 
to time n derives from the fact that each deep essential return originates a bound period 
of length at least ^0 (see Section ET2"j) . Since during the bound periods there cannot be 
any essential return, the number of deep essential returns occurring in a period of length 
n is at most 

2 W 

Proposition 5.2. Given the integers < s < ^ and s < v < n, consider s integers 
7i, . . . , 7 a , each greater than or equal to 0. If Q is large enough, then 

M-^...,») < Qia* 6flj>J. 

Proof. Fix n £ N and take uo £ Vq. Note that the functions v n , s n and & n are constant 
on each ui £ V n . Let u £ uo PI V n be such that v n (u) = v. Then, there is a sequence 
1 < t± < . . . < t v < n of essential return situations. Let Ui denote the element of the 
partition 7\ that contains u. We have ujq D uj\ D . . . D u v = u. For each j £ {0, . . . , v} 
we define the set: 

Qj = U ^' 

and its partition 

Qj = {ujj : UJ £ V n H u) }- 

Let u £ V n be such that & n (u) = s. Then, we may consider 1 < r± < . . . < r s < v with 
indicating that the i-th deep essential return occurs in the r^-th essential return situation. 
Now, set V(0) = Qo = Mo- Fix s integers 1 < r\ < . . . < r s < v. Next, for each j < v we 
define recursively the sets V(j). Although the set V(v) will depend on the fixed integers 
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1 < ?T < . . . < r, < t>, we do not indicate this so that the notation is not overloaded. 
Suppose that V(j — 1) is already defined and rj_i < j < r^. Then, we set 

V{j) = |J u> n - 1/ ) n v(j - 1). 



If j = Ti then we define 



V(j) = |J w n /-^(/ 7i u /_ 7 J n V(j - r 



oj6Q 7 - 



Observe that for every j £ {!,...,«} we have fyT^jji 



< 1. Therefore, we concentrate in 



finding a better estimate for 



Consider that u r . £ Q ri R V(ri) and let cj^-i £ 



<2 ri _i Pi V^(rj — 1) contain u n . We have to consider two situations depending on whether 
t n -\ is an escaping situation or an essential return. 

Let us suppose first that t Ti -\ was an essential return with return depth rj. Then, 



\0J r 



\0J r 



< 

l^n-il |^n-l| 



where cD. r ._i = D f a n {Ui) 



< C 



< C 



fa 1 ifin-l) 

2e" 7 



by m 



-5/3?7 : 



by (E3D- 



Note that when rj_i = r t — 1 then r] = 7i_i. If, on the other hand, rj_i < — 1 then 
is an essential return with depth 77 < 9 < 7i_i. Thus, in both situations, we have 



\0J r 



< 2C- 



-7i 



When i rj _i is an escape situation, instead of using (I3.4|) . we can use (13. 5p and obtain 



< 2C- 



~/3A 



< 2C- 



-5071- 



Observe also that if uJ n -\ ^ 0J n -\ then, because we are assuming that u Ti ^ 0, we have 



> e" 1 - e _e > e" 5 ^- 1 , for large 9. 
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At this point we may write 

lav 



\V{n)\= Y, 



\U ri -l\ 



\U ri -l\ 



w ri eS n nV(rj) 
< 2Ce-^e 5/3 ^- 1 \V{r i - 1)|. 

This yields 

| V(v) | < (2C) S Exp | -(1 - 5/3) Y 7i | e 5/37 ° | V(0) | , 

where 70 is given by the interval uj G Vq. If uj = I( Vo ,k ) with \rj \ > A and 1 < k < t]q, 
then 7 = |r/ |. If uj = (5, 1] or u = [—1, —5), then we can take 70 = 0. 

Now, we have to take into account the number of possibilities of having the occurrence 
of the event V(v) implying the occurrence of the event 7s (^)- The number of possible 
configurations related with the different values that the integers 77, . . .r s can take is 
Hence, it follows that 

M^£...,7>)) <(2cr(^Exp(-(i-5/3)^ 7 4 £ e5/5|70l W 

^ ' I i=l J WoGPo 

<(2C)-QExp j -(1-5/3) J>| (2(1-5)+ £ 

I i=l J \ |r?ol>A 



e 5/3r?o e — 1»? | 



< 3(2C) s Qex P |-(1 - 5/3) , for A large enough 

<QExp|-(l-6^)^ 74 J. 

The last inequality results from the fact that sQ < 2~2i=i 7« an d the freedom to choose a 
sufficiently large 6. □ 

Given the integers < s < s < v < n and the integers 70, 7i, both greater than or 
equal to 0, we consider the event: 



B^ m (n) = < x £ i 7o : v n (x) = v, & n (x) = s and n is a free deep return 

with depth 71 
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Corollary 5.3. Consider the integers 0<s<-^,s<v<n and 70, 71 > 6. If® is large 
enough, then 

* < ( V \exp{-{\ - 6/3)( 7 o + 7i)} • 



The proof of this statement follows easily from Proposition 15. 2\ by observing that, al- 
though n may be an inessential deep return time (instead of an essential deep return), the 
estimates still prevail and for 6 > A large enough we have X] 7 >e e _( ^ 1_6/3 ^ 7 < 1. 

6. The condition D'(u n ) 

Assume that X , X±, . . . is the stationary stochastic process defined in (11. ip with common 
d.f. G a . For z G [—1, 1] the event {Xj(z) > u n } corresponds to the set ((u n , 1]). If n is 
sufficiently large, then /^(wn, 1] = (— a/(1 — u n )/a, — u n )/a) C £/a- We may define 



6 = 0(n) 



1 1 - M r 

-log 

2 a 



(6.1) 



This way, if an exceedance occurs at time j then a deep return with depth over the threshold 
6 must have happened at time j — 1, i.e., if Xj(z) > u n then Xj-x(z) = f{~ l (z) G Uq. 

Remember that the sequence u n is such that n(l — G a (u n )) — > r, as n — > 00, which we 
rewrite as 1 — G a {u n ) x 1/n. Then, by (14. ip . we get u n x 1 — 1/n 2 , which according to 
flBTT]) leads to 

Gxlogn, (6.2) 
meaning that > c, for some c > 0, as n — > 00. 

Observe that we are dealing with very small perturbations of fi for which /|(0) = — 1 
for every j > 2. Thus, one expects that after a deep return to the critical region (a tight 
vicinity of 0) it should take a considerable amount of time before another deep return 
should occur. Since exceedances are related with the occurrence of deep returns then one 
may have a fair amount of belief that condition (12. 2p holds for the sequence Xq, Xi, . . . 

Remark 6.1. If the sequence X ,Xx, . . . was independent, then (12.21) would follow easily 
since 

[n/k] [n/k] [n/k] 

n P{X Q > u n and Xj > u n } — n ^ P{X > u n 
j=i j=i j=i 

< — (1 - G a {u n )) 2 >- ^0. 

K n— >co K k— >oo 

Let us give some insight into the argument we use to prove that (12. 2p holds for Xq, Xi, . . . 

(1) We use the exponential decay of correlations (see (13.61) ) to compute a turning instant 
T = T(n) such that the dependence between X Q and Xj with j > T is negligible. 
This suggests the splitting of the time interval {1, . . . , [71/ A;]} into {1, . . . , T} and 
{T + 1, . . . ,[n/k]}, when n is sufficiently large. 
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(2) During the time interval {T + 1, . . . ,[n/k]} we use the fact that for j > T the r.v. 
Xj is almost independent of Xq and argue like in Remark 16.11 

(3) For j G {1, . . . , T} we appeal to Corollary 15.31 to bound P{X > u n and Xj > u n } 
and afterwards we use the fact that, for n large, T [n/k] to finish the proof. 

Step (1) 

Taking <p = ijj = l( Un ,i] in (JUSj), we get 

\u a {X > u n and Xj > u n } - [fi a {X > u n }] 2 \ = 

• l(«„,i] ° Pa.dn a - I l (Wre! i]d/i a 



where we may assume that C is the same for all n G N, because ||l( Un ,i]||oc = 1 an d the 
total variation of l( u „,ii is equal to 1, for every n G N. 

We compute T = T(n) such that for every j > T we have 

1 



CV < 



n 



3 ' 



Since CV < ^3 <^ j > loo .^-i (3 logn + logC), we simply take, for n sufficiently large, 



r 



log c 



— logn. 



(6.3) 



For fixed k and n sufficiently large, we have that T < [n/k]. Hence, we may write 



[n/k] 



n 



P{X > u n and Xj > u n } 



3=1 



[n/k] 



= n P{X > u n and Xj > u n } + n P{X > u n and Xj > u n }. 

j=\ j=T+l 

In step (2) below we deal with the second term in the sum, leaving the first term for step 
(3). 



Step (2) 



Let us show that lim sup^^ n $^=t+i P{Xq > u n and Xj > u n } — > 0, as k — > 00. By 
choice of T we have 



i n / k ] / j\ 

^ P{X > u n and Xj > n n } < n I (1 - G a (n n )) 2 + — 1 [n/k] 



j=T+l 
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Now, £(1 - G a {u n )f + & ► £ ► and the result follows. 

Step (3) 

We are left with the burden of controlling the term n X)J=i P{Xq > u n and Xj > «„}. 
We begin with the following lemma that will enable us to bound the number of exceedances 
occurring during the time period {1, . . . , T}. In what follows, we are always assuming that 
n is large enough so that G > A. 

Lemma 6.2. If a deep return occurs at time t (with depth over the threshold Q), then the 
next deep return can only occur after t + 0/2. 

Proof. For every x G Uq, the bound period associated to x is such that p{x) > 0/2, by 
Section EJ ©• For all j < [6/2] we have 

Am 

\fl(x)\ > \fi(0) \ - e~ pj > e~ a ^ - e~? j > e~ aj (l - e^"^') 

> e~ aj (1 - e (a_/3) ) , since a - < 

> e - Q0/2 (1 - e^"^) , since j < 9/2 

> e~ o0 , if n is large enough so that 1 — e Q_ ^ > e~ a& ^ 2 

> e _e , since a < 1. 

□ 

As a consequence of Lemma 16.21 we have that the maximum number of exceedances up 
to time T is at most 2T/0. 

Next lemma shows that in the bound period following a deep return over the level 0, 
there cannot occur deep bound returns before time T. 

Lemma 6.3. For every a e BC , if n is sufficiently large then for every x 6 Uq we have 
that fl(x) (jz Uq for all j < min{p(x),T}. 

Proof. Consider the map h : [l,oo) —> R given by h(y) = e~ a Vv — e _/3y . There is J = 
J(a, 0) e N such that ti{y) < for every y > J. Let D = e""^. 

Having in mind that by f)6.2p and ( 16. 3 p we have = 0(n) x logra and T x logn, 
respectively, assume that n is large enough so that: 

(1) e- < D/2- 

(2) by continuity of f a , for every x G Uq and all j < J we have \fl(x) — /„(0)| < D/2; 

(3) e~ aVT - e-^ T > e" e . 

Let x G Uq and consider j < min{p(a;), J}. By assumption (j2J) we have \fl(x) — fi(0)\ < 
D/2. Using fIBAp and hypothesis §T§ we get 

\fi(x)\ > \fim - \fi(x) - fim > d-d/2 > d/2 > e ^ e . 
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Now, let x G Uq and consider J < j < min{p(x), T}. Using the definition of bound 
period, (IBAj) and condition ([3]) we have 



\fi(x)\>\fi(o)\-\fi(x)-fi 



where the third inequality derives from the fact that on [J, oo) the function h is decreasing. 

□ 

As a consequence, if there is a deep return at time t, then we cannot have bound 
returns during the time period (t, t + T] . This way, we may use Corollary 15.31 to estimate 

\u @ nf-m @ \JoTj<T. 

Observe that for n large enough, we have 

P{X > u n and Xj > u n } = \i a {(u n , 1] n f~ J (u n , 1]} 

= A*.{/a~ 1 («n,l]n/- &+1) («n,l]}, 

and fa X {Uni 1] n /~ J ' {f~ x {u n , 1]) C Uq PI f~iJJ®. In what follows we will use "const" to 
denote several positive constants independent of n. 
Note that 

uenf^uec u U U U 

s=0 ti=S7o>©7i>0 



Hence, by Corollary 15.31 we have 

2j/e j , x 2T/e t 



% n /-f/el <EEEE (")e-"-«")«— » < const £ £ (') 

s=0 «=s 7 >9 7i>e ^ ' s=0 v=° ^ ' 

2T/B T , , 2T/6 , . 

< const ]T ^ I V^ 20 < const £ tI )e 



V ) e -(l-6/3)26 



At this point, we estimate 2T/0. Recalling (16.21) and (16. 3p . one easily gets that there 
exists a constant Ci > such that 2T/6 < Ci, for n sufficiently large. So, to proceed with 
the estimation \Uq D J~^Uq\, assume that n is sufficiently large so that 2T/6 < Ci and 
d < T. Then, 

2T/0 . N 2T/0 



e -(l-6 / 8)2G 



|C/© n /-^ e | < const £ t( V^ 6 « 20 < const £ j 

s=0 \ S ' s=0 ^ X ' 

< const^T^e"^ 20 < constT Cl+1 e-( 1 - 6 ^ 20 . 

So far we managed to obtain an estimate for \Uq D /~ J 'Z7e|. Let us see how too derive 
one for P{X > u n and Xj > u n } < \i a {Uq D /^Uq}. By assertion ([3]) in Section we 
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have 

Ma {Uq n f- j U e } = [ P 1{x)dx + / P a 2 (x)dx 

Ju e nf- J u e Ju e nf- J u @ 



<const.\U e nf- j Ue\+y2 I 

.7=1 Ju 



e 2 



Using (jES} we get |/-['(0)| > e" ae , for all j < 6 2 . Thus, for all x G U and 9 



large enough, we have y\x — fl(0)\ > y/e Q0 — e > 1/2. e qQ / 2 ; which implies that 

-, 1 . < 2e a0 / 2 < 2e /3e (recall that /3 > a). Consequently, for sufficiently large 6, 

V \x-fi(o)\ 

{U e n /~ i E/'e} < const. \U e n /~ i C/'e| + const. |£/© n /~ J 't/'e|.e^ e + const.e° 2logL9 
< const.T Cl+1 e-( 1 - 7/3 ) 20 . 

Now, since e -0 < const (1 — G a (u n )), we finally conclude: 

T T 

n P{X > u n and Xj > u n } < const • n T Cl+1 (l - G a {u n )) 2{l - w) 
j=i j=i 

< const • nT Cl+2 (l - G a {u n )) 2 ^ 7 ^ 

< const • n(logn) Cl+2 (l - G a (n n )) 2 " 14/3 

< const • [n(l - G a (n n ,))] 3 / 2 (1 - G a (n n )) 1/2 ~ 14/3 , 
for sufficiently large n. The result follows because 

hm [n(l - G a K))] 3 / 2 (l - CW) 1 ^ = r 3 / 2 -0 = 0. 

7. Simulation Study 

In this section we present a small simulation study illustrating the finite sample behavior 
of the normalized M n , defined in (ll.2|) . for the Benedicks-Carleson quadratic map /^(a?) = 
1 — 2x 2 . We have that P{Xo < x} = (^(x) = l/2 + arcsin(x)/7r. According to Theorem IA1 
the normalizing sequence is a n — (1 — cos(7r/n)) _1 , for each n G N, and the theoretical 
limiting distribution for a n {M n — 1) is 



e-^ 172 ,x<0 
1 ,x > 



We performed the following experiment. We picked at random (according to the d.f. 
Gq) a point z in the interval [—1,1], computed its orbit up to time n and calculated 
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M n (z) = maxjz, f 2 (z), . . . , f% (z)}- We repeated the process m times to obtain a sample 
{M n (zi), . . . , M n (z m )} and approximated, for certain values of x, P{a n (M n — b n ) < x} by 



^ m 

i=l 



(zi)-b n )<x}i 



(7.1) 



, - jl if a»(M n (zi) - 6 n ) < x . 

where l {an{Mn{Zt) _ bn) < x} = < , for each 1 < i < m. 

\0 if a n (M n (zi) - bn) > x 

In Table [T]we present the results obtained by realizing the above experiment, considering 
different values of x, n and m and compare them with the theoretical limiting ones given 
by H(x). 







n = 


1000 


n = 10000 


n = 20000 


X 


H(x) 


m = 1000 


m = 10000 


m = 10000 


m = 20000 


m = 20000 


-0.001 


0.9689 


0.976 


0.9671 


0.9719 


0.9677 


0.9708 


-0.01 


0.9048 


0.894 


0.9079 


0.9056 


0.9076 


0.9052 


-0.1 


0.7289 


0.724 


0.7263 


0.7323 


0.7303 


0.7335 


-0.3 


0.5783 


0.569 


0.5773 


0.5823 


0.5840 


0.5813 


-0.5 


0.4931 


0.491 


0.4906 


0.5012 


0.4984 


0.4941 


-0.7 


0.4332 


0.407 


0.4272 


0.4403 


0.4374 


0.4338 


-1 


0.3679 


0.388 


0.3631 


0.3663 


0.3729 


0.3678 


-3 


0.1769 


0.164 


0.1731 


0.1748 


0.1833 


0.1729 


-5 


0.1069 


0.124 


0.1024 


0.1092 


0.1108 


0.1056 


-8 


0.0591 


0.059 


0.0510 


0.0557 


0.0617 


0.0580 


-10 


0.0423 


0.049 


0.0350 


0.0435 


0.0438 


0.0414 


-30 


0.0042 


0.002 


0.0031 


0.0033 


0.0048 


0.0041 


-50 


0.00085 


0.001 


0.0007 


0.0009 


0.0007 


0.0009 


Table 1. Simul 


ation results 



As one may verify the results of the experiment are quite close to the asymptotic theoret- 
ical ones and there is a general tendency of getting better as n increases which is precisely 
the behavior we were expecting. It is also noticeable that there is an improvement when m 
increases, which is also predictable since our approximation (17.11) gets to be more accurate. 
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